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Abstract 

Cosmological shock waves are induced during hierar- 
chical formation of large-scale structure in the universe. 
Like most astrophysical shocks, they are collisionless, 
since they form in the tenuous intergalactic medium 
through electromagnetic viscosities. The gravitational 
energy released during structure formation is trans- 
ferred by these shocks to the intergalactic gas as heat, 
cosmic-rays, turbulence, and magnetic fields. Here we 
briefly describe the properties and consequences of the 
shock waves in the context of the large-scale structure 
of the universe. 

Keywords Cosmic-rays • Large-scale structure of uni- 
verse ■ Magnetic fields • Shock waves • Turbulence 



1 Introduction 

Shock waves are ubiquitous in astrophysical environ- 
ments; from solar winds to the largest scale of the 



universe ( Rvu et al. 20031 ). In the current paradigm 



of the cold dark matter (CDM) cosmology, the large- 
scale structure of the universe forms through hierarchi- 
cal clustering of matter. Deepening of gravitational po- 
tential wells causes gas to move supersonically. Cosmo- 
logical shocks form when the gas accretes onto clusters, 
filaments, and sheets, or as a consequence of the chaotic 
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flow motions of the gas inside the nonlinear structures. 
The gravitational energy released during the formation 
of large-scale structure in the universe is transferred by 
these shocks to the intergalactic medium (IGM). 

Cosmological shocks are collisionless shocks which 
form in a tenuous plasma via collective electromagnetic 
interactions between bar yonic particles and electromag- 
netic fields (|QuestlE"988h . They play key roles in gov- 
erning the nature of the IGM through the following pro- 
cesses: in addition to the entropy generation, cosmic- 
rays (CRs) are produced via diffusive shock accelera - 
tion (DSA) |Bel3ll978t iBlandford and Ostrikerlll978h . 
magnetic fie lds are generated via the Biermann battery 
mechanism ( Kulsru d et al. I ll997h and Weibel nstability 
dMedvedey et al . 2006) and also amplified by streaming 
CRs (IBelll 12004 ) , and vorticity is generated at curved 
shocks (|Binnevlll974h . 

Cosmological shocks in the intergalactic space have 
been studied in details using various hydrodynamic sim- 
ulations for the cold dark m atter cosmology with cos- 

molqgical constant (ACD M) ( Ryu et al~ll2003 ; Pfrommer et al.^ 
2006t iKang et al.l l2008h . In this contribution, we de- 



scribe the properties of cosmological shocks and their 
implications for the intergalactic plasma from a sim- 
ulation us ing a PM/Eule rian hydrodynamic cosmol- 
ogy code ( Rvu et al. 1993) with the following param- 
eters: n BM = 0.043, S1 DM = 0.227, and Q A = 0.73, 
h = iV(!00 km/s/Mpc) = 0.7, and a 8 = 0.8. A 
cubic region of comoving size 100 h~ x Mpc was sim- 
ulated with 1024 3 grid zones for gas and gravity and 
512 3 particles for dark matter, allowing a uniform spa- 
tial resolution of Al = 97.7/i _1 kpc. The simulation 
is adiabatic in the sense that it does not include ra- 
diative cooling, galaxy/star formation, feedbacks from 
galaxies/stars, and reionization of the IGM. A temper- 
ature floor was set to be the temperature of cosmic 
background radiation. 
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Fig. 1 Two-dimensional images showing x-ray emissivity (top left), locations of shocks with color-coded shock speed V s 
(top right), perpendicular component of vorticity (bottom left), and magnitude of vorticity (bottom right) in the region of 
(25 /i _1 Mpc) 2 around a galaxy cluster at present (z = 0). Color codes V a from 15 (green) to 1,800 km s" 1 (red). 



2 Properties of Cosmological Shocks 

As a post-processing step, shocks in the simulated vol- 
ume are identified by a set of criteria based on the shock 
jump conditions. Then the locations and properties of 
the shocks such as shock speed (V s ), Mach number (M), 
and kinetic energy flux (/kin) are calculated. 

In the top panels of Figure 1, we compare the loca- 
tions of cosmological shocks with the x-ray emissivity 
in the region around a cluster of galaxies, both of which 
are calculated from the simulation data at redshift 
z = 0. External shocks encompass this complex nonlin- 
ear structure and define the outermost boundaries up 
to ~ 10 hr 1 Mpc from the cluster core, far beyond the 
region observable with x-ray of size ~ 1 hT 1 Mpc. In- 



ternal shocks are found within the region bounded by 
external shocks. External shocks have high Mach num- 
bers of up to M ~ 10 3 due to the low temperature of 
the accreting gas in the void region. Internal shocks, 
on the other hand, have mainly low Mach numbers of 
M < 3, because the gas inside nonlinear structures has 
been previously heated by shocks and so has high tem- 
perature. 

The frequency of cosmological shocks in the simu- 
lated volume is represented by the quantity 5, the area 
of shock surfaces per unit comoving volume, in other 
words, the reciprocal of the mean comoving distance 
between shock surfaces. In the top left panel of Figure 
2, we show S(V S ) per unit logarithmic shock speed in- 
terval at z = 0. We note that the frequency of low speed 



Shock Waves in the Large-Scale Structure of the Universe 



3 



„ -l 

£ -2 - 

§ ; 

-3 L 

to : 




0.6 



12 3 
log V, (km/s) 




1.5 2 2.5 3 3.5 
log V s (km/s) 




1.5 2 3. 5 3 3.5 
log V s (km/s) 



1.5 a 3. 5 3 3.5 
log V s (km/s) 



Fig. 2 (Top left) Reciprocal of the mean comoving distance 
between shock surfaces at z — in units of l/(fa _1 Mpc). 
(Top right) Kinetic energy flux passing through shock sur- 
faces per unit comoving volume at z = in units of 10 40 ergs 
s _1 (/i _1 Mpc)~ 3 . (Bottom left) Thermal energy dissipated 
(dotted line) and CR energy generated (solid line) at shock 
surfaces, integrated from z = 5 to 0. (Bottom right) Cu- 
mulative energy distributions. The energies in the bottom 
panels are normalized to the total thermal energy at z — 0. 
All quantities are plotted as a function of shock speed V a . 



shocks with V s < 15 km s is overestimated here, 
since the temperature of the intergalactic medium is 
unrealistically low in the adiabatic simulation without 
the cosmological reionization process. Although shocks 
with V s ~ a few x 10 km s _1 are most common, those 
with speed up to several xlO 3 km s _1 are present at 
z = 0. The mean comoving distance between shock sur- 
faces is 1/S ~ 3 /i -1 Mpc when averaged over the entire 
universe, while it is ~ 1 /i Mpc inside the nonlinear 
structures of clusters, filaments, and sheets. 

In order to evaluate the energetics of cosmological 
shocks, the incident shock kinetic energy flux, /km = 
(l/2)piV s 3 , is calculated. Here pi is the preshock gas 
density. Then the average kinetic energy flux through 
shock surfaces per unit comoving volume, F, is calcu- 
lated. The top right panel of Figure 2 shows F(V S ) 
per unit logarithmic shock speed interval. Energeti- 
cally the shocks with V s > 10 3 km s _1 , which form in 
the deepest gravitation potential wells in and around 
clusters of galaxies, are most important. Those respon- 
sible for most of shock energetics are the internal shocks 
with relatively low Mach number of M ~ 2 — 4 in the 
hot IGM, because they form in the high-density gas in- 
side nonlinear structures. On the other hand, external 
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Fig. 3 Gas thermalization efficiency, S(M), and CR ac- 
celeration efficiency, r/(M), as a function of Mach number. 
Symbols are the values estimated from numerical simula- 
tions based on a DSA model and dotted and dashed lines 
are the fits. Solid line is for the gas thermalization efficiency 
for shocks without CRs. 



shocks typically form in accretion flows with the low- 
density gas in voids, so the amount of the kinetic energy 
passed through the external shocks is rather small. 



3 Energy Dissipation at Cosmological Shocks 

In addition to the gas entropy generation, the accelera- 
tion of CRs is an integral part of collisionless shocks, in 
which electromagnetic interactions between plasma and 
magnetic fields provide the necessary viscosities. Supra- 
thermal particles are extract ed from the shock-heate d 
thermal particle distribution ( Malkov and Drurvll2001 ). 
With 10~ 4 — 10~ 3 of the particle flux passing through 
the shocks injected into the CR population, up to ^60% 
of the kinetic energy of strong quasi-parallel shocks can 
be converted into CR ions and the nonlinear feedback to 



the u nderlying flow can be substantial (jKang and Jones 
2005). At perpendicular shocks, however, the CR in- 
jection and acceleration are expected to be much less 
efficient, compared to parallel shocks, since the trans- 
port of low energy particles normal to the average field 
direction is suppressed. So the CR acceleration depends 
sensitively on the mean magnetic field orientation. 

Time-dependent simulations of DSA at quasi-parallel 
shocks with a thermal leakage injection model and a 
Bohm-typc diffusion coefficient have shown that the 
evolution of CR modified shocks becomes self-similar, 
after the particles are accelerated to relativistic en- 
ergies and the pr ecursor compression re a ches a time- 



asymptotic state (jKang and Joned I2005L 120071 ). The 
self-similar evolution of CR modified shocks depends 
somewhat weakly on the details of various particle- wave 
interactions, but it is mainly determined by the shock 
Mach number. Based on this self-similar evolution, we 
can estimate the gas thermalization efficiency, S(M), 
and the CR acceleration efficiency, rj(M), as a function 



1 



shock Mach number, which represent the fractions of 
the shock kinetic energy transferred into the thermal 
and CR energies, respectively. Figure 3 shows the re- 
sults of such DSA simulations (jKang et al.ll2008f h From 
the figure, we expect that at weak shocks with M < 3 
the energy transfer to CRs should be < 10 % of the 
shock kinetic energy at each shock passage, whereas at 
strong shocks with M > 3 the transfer is very efficient 
and the flow should be significantly modified by the CR 
pressure. 

Adopting the efficiencies shown in Figure 3, the ther- 
mal and CR energy fluxes dissipated at each cosmolog- 
ical shock can be estimated as S(M)fki n and T}{M)fun, 
respectively. Then the total energies dissipated through 
the surfaces of cosmological shocks during the large- 
scale structure formation of the universe can be cal- 
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culated (|Rvu et al.ll2003t iKang et all 120081 ). The bot- 
tom panels of Figure 2 show the thermal and CR en- 
ergies, integrated from from z — 5 to 0, normalized 
to the total gas thermal energy at z — 0, as a func- 
tion of shock speed. Here we assume CRs are freshly 
injected at shocks without a pre-existing population. 
The shocks with V s > 10 3 km s" 1 are most responsible 
for the shock dissipation into heat and CRs. The fig- 
ure shows that the shock dissipation can count most o f 
the gas thermal energy in the IGM (jKang et alJl2005 ). 
The ratio of the total CR energy, Y"cr(> V s , m i n ), to the 
total gas thermal energy, Y t h{> V,, ,„;„), dissipated at 
cosmological shocks throughout the history of the uni- 
verse is about 0.4 (where V Stia \ n is the minimum shock 
speed), giving a rough estimate for the energy density 
of CR protons relative to that of thermal gas in the 
IGM as £cr p ~ 0.4e t h crm • 

Because of uncertainties in the DSA model such as 
the field obliquity and the injection efficiency, however, 
it is not meant to be an accurate estimate of the CR 
energy in the IGM. Yet the results imply that the IGM 
could contain a dynamically significant CR population. 



4 Turbulence Induced by Cosmological Shocks 

Vorticity can be generated in the IGM either directly 
at curved cosmological shocks or by the baroclinity of 
flow. The baroclinity is resulted mostly from the en- 
tropy variation induced at cosmological shocks. There- 
fore, the baroclinic vorticity generation also can be 
attributed to the presence of cosmological shocks. A 
quantitative estimation of the vorticity in the IGM was 
made using the data of the simulation described in §1 
(|Rvu et al.ll2008l ). The bottom panels of Figure 1 show 
the distribution of the vorticity around a cluster com- 
plex. The distribution closely matches that of shocks, 
as expected. 
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Fig. 4 Volume fraction with given temperature and vortic- 
ity magnitude (top left), with given gas density and vorticity 
magnitude (top right), with given temperature and mag- 
netic field strength (bottom left), and with given gas den- 
sity and magnetic field strength (bottom right) at present 
[z = 0). Here, t agc is the present age of the universe. 



The top panels of Figure 4 show the magnitude of 
the vorticity in the simulated volume as a function of 
gas temperature and density. Here the vorticity mag- 
nitude is given in units of the reciprocal of the age of 
the universe, l/t age - There is a clear trend that the 
vorticity is larger in hotter and denser regions. At the 
present epoch, the rms vorticity is ^rms^agc ~ 10 to 30 
in the regions associated with clusters/groups of galax- 
ies (T > 10 7 K) and filaments (10 5 < T < 10 7 K), 
whereas it is on the order of unity in sheetlike struc- 
tures with 10 4 < T < 10 5 K and even smaller in voids 
with T < 10 4 K. 

With Wdy = 1 /<■<-> interpreted as the local eddy 
turnover time, u> x i age represents the number of eddy 
turnovers of vorticity in the age of the universe. It takes 
a few turnovers for vorticity to decay into smaller ed- 
dies and develop into turbulence. So with Lo rms t agc ~ 
10 — 30, the flows in clusters/groups and filaments is 
likely to be in the state of turbulence. On the other 
hand, with w rms t agc ^ 1 the flow in sheetlike structures 
and voids is expected to be mostly non-turbulent. 

In order to estimate the energy associated with the 
turbulence in the IGM, the curl component of flow mo- 
tions, Ucurii which satisfies the relation V x v cul -i = 
V x v, is extracted from the velocity field. As vor- 
ticity cascades and develops into turbulence, the en- 
ergy (l/2)/5tr , is transferred to turbulent motions, 
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so it can be regarded a s the turbulence energy, eturb 



As shown in iRvu et al.l (|2008l ) , eturb < etherm in clus- 



ters/groups. In particular, the mass-averaged value 
is (eturb/ethcrm) mass = 0.1 - 0.3 in the intracluster 
medium (ICM), which is in good agreement with the 
observationally inferre d values in cluster core regions 
(Schucckcr ct al . 2004). '. n filaments and sheets, this 
ratio is estimated to be 0.5 < (eturb/etherm) < 2 and it 
increases with decreasing temperature. 



5 Intergalactic Magnetic Field 

How have the intergalactic magnetic fields (IGMFs) 
arisen? The general consensus is that there was no 
viable mechanism to produce strong, coherent mag- 
netic fields in the IGM prior to the formation of 
large -scale structure and galaxies ( Kulsrud and Zweibell 
2008). However, it is reasonable to assume that weak 
seed fields were created in the early universe. A num- 
ber of mec hanisms, including t he Biermann battery 
mec hanism (IKulsrud et alJll997l ) and Weibel instabil- 
ity ( Medvedev et al. 20061 ) worki ng at early cosmologi- 
cal sh ocks, have been suggested ( Kulsrud and Zweibell 
2008). The turbulence described in §4 then can am- 
plify the seed fields in the IGM through the stretch- 
ing of field lines, a process known as the turbulence 
dynamo. In this scenario the evolution of the IGMFs 
should go through three stages: (i) the initial exponen- 
tial growth stage, when the back-reaction of magnetic 
fields is negligible; (ii) the linear growth stage, when the 
back-reactio n starts to operate; and (iii) the final satu - 
ration stage l|Cho and Vishniadl200oHcho et al.ll2008h . 

In order to estimate the strength of the IGMFs re- 
sulted from the dynamo action of turbulence in the 
IGM, we model the growth and saturation of magnetic 
energy as: 



SB 



0.04 x cxp[(<' - 4)/0.36] 



for t' < 4, 



-turb 



= { (0.36/41) x (f - 4) + 0.04 for 4 < t' < 45, 



0.4 



for t' > 45, 



based on a simulatio n of incompressible magnetqhydro - 
dynamic turbulence ( Rvu et al. 2008 : Cho et al. 2008). 
Here t' = i/i e ddy is the number of eddy turnovers. This 
provides a functional fit for the fraction of the turbu- 
lence energy, Eturb, transfered to the magnetic energy, 
Sbi as a result of the turbulence dynamo. 

The above formula is convoluted to the data of the 
simulation described in §1, setting t' = u> x t age , and the 
strength of the IGMFs is calculated as B = (8-keb) 1 ^ 2 - 
The resulting magnetic field strength is presented in 
the bottom panels of Figure 4. On average the IGMFs 
are stronger in hotter and denser regions. The strength 



of the IGMFs is B > 1/iG inside clusters/groups (the 
mass-averag ed value for T > 10 7 K), ~ 0.1/xG around 
clusters/groups (the volume-averaged value for T > 10 7 
K), and - 10 nG in filaments (with 10 5 < T < 10 7 
K) at present. The IGMFs should be much weaker 
in sheetlike structures and voids. But as noted above, 
turbulence has not developed fully in such low density 
regions, so our model is not adequate to predict the 
field strength there. 

We note that in addition to the turbulence dynamo, 
other processes such as galactic winds driven by super- 
nova explosions and jets from active galactic nuclei can 
further streng then the magnet ic fields to the IGM (for 
references, see lRvu et al.l l2008) 



6 Conclusion 

Shocks are inevitable consequences of the formation 
of large-scale structure in the universe. They heat 
gas, accelerate cosmic-ray particles, produce vortic- 
ity and turbulence, and generate and amplify mag- 
netic fields in the IGM. By applying detailed mod- 
els of the DSA and the turbulence dynamo to the 
data of a cosmological hydrodynamic simulation of a 
concordance ACDM universe, we have made the fol- 
lowing quantitative estimates: ECRp/stherm ~ 0.4 in 

the IGM, <£turb/£the r m) m as S = 0.1 - 0.3 ill the ICM, 

(Eturb/^therm) ~ 0.5 — 2 in filaments and sheets, B > 
1/iG inside clusters/groups, B ~ 0.1/iG around clus- 
ters/groups, and B ~ 10 nG in filaments at present. 
Our results suggest that the non-thermal components 
can be energetically significant in the intergalactic 
plasma of large-scale structure, as in the interstellar 
plasma inside Our Galaxy. 
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